Kardar-Parisi-Zhang interfaces bounded by long-ranged potentials 
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We study unbinding transitions of a non-equilibrium Kardar-Parisi-Zhang (KPZ) interface in the 
presence of long-ranged substrates. Both attractive and repulsive substrates, as well as positive and 
negative Kardar-Parisi-Zhang nonlinearities, are considered, leading to four different physical situ- 
ations. A detailed comparison with equilibrium wetting transitions as well as with non-equilibrium 
unbinding transitions in systems with short-ranged forces is presented, yielding a comprehensive 
picture of unbinding transitions, and of their classification into universality classes. These non- 
equilibrium transitions may play a crucial role in the dynamics of wetting or growth of systems with 
intrinsic anisotropics. 
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I. INTRODUCTION 

Spatial constraints in systems where two (or more) 
bulk phases coexist may lead to wetting transitions. This 
is the case, for example, of confined fluids where one of 
two coexisting equilibrium phases (the liquid, say) is in 
contact with a substrate with an interface separating it 
from the second phase (the gas) at infinity. The liquid 
does not wet the substrate if the thickness of the liq- 
uid film is finite (there is a microscopic quantity of liq- 
uid). On the other hand, the substrate is wet if there 
is a macroscopically thick liquid film on it. A wetting 
transition is said to occur when the substrate changes 
from not being wet by the liquid to being wet. Typi- 
cally, two types of wetting transitions can be considered: 
by increasing the temperature at bulk coexistence one 
may find either critical wetting or a discontinuous transi- 
tion; by varying the chemical potential while the temper- 
ature is fixed above the wetting transition temperature 
one finds a complete wetting transition, at bulk coexis- 
tence. Under equilibrium conditions, a completely anal- 
ogous transition (often called drying) may occur when 
the substrate preferentially adsorbs the gas phase 1]. 

Effective interfacial potentials are useful coarse-grained 
models that have played a key role in understanding 
a large variety of equilibrium wetting problems P, 01 ■ 
These potentials, V{h), are functionals of the interfa- 
cial height (measured from the substrate), /i(x). In this 
framework, wetting transitions are described as the un- 
binding of the (say liquid-vapor) interface from the sub- 
strate, with the effective binding potential determined 
by the microscopic forces between the constituents of the 
substrate and those of the bulk phases. Typically, ex- 
ponentials and power-law decaying potentials V{h) have 
been considered for systems dominated by short-ranged 
and long-ranged forces, respectively. 

There exist a large amount of phenomena describable 
in terms of equilibrium wetting, either under short-range 
or under long-ranged interactions, while it has only re- 



cently been recognized that non-equilibrium effects, such 
as anisotropies in the interface growing rules, may play a 
crucial role in describing some experimental situations. 
Within this perspective, short-ranged non- equilibrium 
wetting has been studied and some interesting novel 

phenomenology has been elucidated (see for recent re- 
views). In particular, liquid-crystals molecular-beam 
epitaxial systems, as GaAS |3, or materials exhibiting 
Stranski-Krastanov instabilities 01, appear to be good 
candidates to require a non-equilibrium wetting descrip- 
tion. However, some of these systems, as well as many 
others not enumerated, might include effective long-range 
substrate-interface effects as also occurs in the equilib- 
rium case. 

Our goal in this paper is to fill this gap by providing a 
general and systematic theory of non-equilibrium wetting 
under the presence of effective long-ranged interactions. 
First, we briefly review the equilibrium situation to set 
up the theoretical framework and, afterwards, generalize 
it to embrace non-equilibrium situations. 

In equilibrium, two types of analytical approaches are 
available: static studies based on the ensemble theory 
Q and dynamical, stochastic approaches that allow in- 
vestigating relaxational aspects. The second approach is 
amenable to non-equilibrium extensions and is the one we 
employ. Thus, consider the simple Edwards- Wilkinson 
dynamics 9] subject to a bounding force (i.e. the deriva- 
tive of the bounding potential) lOj: 



dth{yi,t) = V'^h 



dVjh) 
dh 



+ ari{^,t). (1) 



This includes (i) the usual diffusion term, computed as 
minus the derivative of a standard surface-tension term, 
(ii) the driving force, a, related to the chemical poten- 
tial difference between the two phases, (iii) the Gaussian 
white noise, ?7(x, i), and (iv) the bounding force, which 
may derive from a short-ranged potential 



V{h) = 



-ph 



-qh 



(2) 



2 



or from a long-ranged one 



^ ' phP qhi 



(3) 



where, 6, c > 0, and p < q are parameters. This last 
form, Eq.Q, is known to be the correct functional form 
for systems where the molecules interact through van der 
Waals forces [Tll |. 

By varying the chemical-potential, a, one controls the 
average interfacial distance from the wall: small for 
a < flc (non-wet phase), large for a « Oc, and increasing 
steadily with time for a > Oc (wet phase), i.e. the system 
exhibits an unbinding transition at a — ac- The inter- 
face potentials V{h) are, in all cases, harshly repulsive 
at small h to model the impenetrability of the substrate. 
The parameter b vanishes linearly with the temperature, 
at the (mean-field) critical wetting temperature, and rep- 
resents the affinity or preference of the substrate for one 
of the bulk phases (usually the liquid) . We consider three 
distinct situations (see Fig. QJ: 

1. Repulsive potential: complete wetting. If 6 > the 
potential describes the presence of a bounding substrate 
alone. In this case, the broken symmetry induced by the 
substrate leads to the divergence of the average position 
of the interface, at coexistence, Oc = 0: i.e. the system 
undergoes a complete wetting transition. 

The latter is described by Eq. Q with the potential 
taken from Eq. Two different regimes depending on 
the value of p have been reported: for p < 2 mean-field 
scaling holds and (h) ^ t^/(p+^)^ while if p > 2 fluc- 
tuations take over and the velocity is controlled by the 
intrinsic roughness of a free Edwards- Wilkinson, lead- 
ing to a fluctuation-dominated regime characterized by 
(h) ~ i^/^ (0, or logarithmic growth for two-dimensional 
interfaces) . These results are derived in a formal way and 
extended in the appendix. 

2. Attractive potential: first order unbinding. For b < 0, 
by contrast to the complete wetting case, the surface does 
not promote the growth of the liquid phase and conse- 
quently there is no wetting phase even at bulk coexis- 
tence, a = 0. V{h) exhibits a local minimum near the 
substrate, that binds the interface in the presence of ther- 
mal fluctuations, and the width of the wetting layer is fi- 
nite (microscopic) at a = 0. We may, however, observe a 
first-order unbinding transition that occurs as a changes 
from positive (stable bulk liquid) to negative (stable bulk 
gas) values. 

3. Critical wetting. At a particular value of 6 = 6c (&c ~ 
in mean-field but more generally be is small and nega- 
tive) critical wetting may be observed, with a character- 
istic non-trivial phenomenology. This situation requires 
the fine tuning of two independent parameters {b — b^, 
a = Oc). This critical transition is more difficult to treat 
theoretically and less likely to be found in real systems 
and will not be discussed in this paper. 
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FIG. 1: Effective potentials as derived from Eq. @ in the 
pinned (a < 0) and the depinned (a > 0) phases, and at co- 
existence (a = 0). (a) repulsive walls {b, c > 0); (b) attractive 
walls (6 < 0, c> 0). 



The best way to extend equilibrium approaches to 
more general, non-equilibrium, situations is to consider 
the simplest and widely studied non-equilibrium exten- 
sion of the Edwards- Wilkinson equation, i.e. the Kardar- 
Parisi-Zhang (KPZ) llSj interfacial dynamics in 
the presence of effective bounding potentials, as the ones 
we have described before. This strategy has been fol- 
lowed in a series of recent papers for systems with short- 
ranged (attractive and purely repulsive) potentials 
and will be extended in the present work to the case 
of long-ranged potentials. We will discuss the phase dia- 
grams for both purely repulsive and attractive potentials, 
paying special attention to criticality and to the compar- 
ison with equilibrium wetting and non-equilibrium short- 
ranged unbinding. We will focus on one-dimensional in- 
terfaces (separating two-dimensional bulk phases), and 
mention briefly two-dimensional interfaces in the conclu- 
sions. 

The paper is organized as follows. In section 2 we intro- 
duce the non-equilibrium unbinding model. In section 3, 
we review known results for non-equilibrium short-ranged 
unbinding. Section 4 contains the main body of the pa- 
per, including both analytical and numerical results for 
purely repulsive and attractive potentials. Finally, the 
main conclusions are presented together with a discus- 
sion of our results. 



II. NON-EQUILIBRIUM LONG-RANGED 
UNBINDING: THE MODEL 

Our model consists in a KPZ non-equilibrium inter- 
face 9, 13J in the presence of a long-ranged, bounding 
potential Eq. 

dth = V^h + \{Vhf +a+j^ + -^+ a77(x, t), (4) 



where A 7^ is the coefficient of the non-linear KPZ term, 
the only new ingredient added to the equilibrium wetting 



3 



Langevin Eq. 

Note that in equihbrium the time-dependent probabil- 
ity distribution P{h, t) is symmetric for the free interface 
and therefore it does not make any difference which side 
faces the substrate. By contrast, under non-equilibrium 
conditions, owing to the — > — /i asymmetry of the KPZ- 
equation, it depends on the sign of A that the substrate 
probes either one tail or the other of a KPZ probabil- 
ity distribution that is no longer symmetric. Thus, for a 
given bounding potential two different situations must be 
considered. Therefore, we will investigate systems with 
positive and negative values of A (without loss of gener- 
ality we take A = ±1), and with both attractive (6 < 0) 
and repulsive {b > 0) potentials, i.e. we consider four 
distinct cases. The focus is mainly on one-dimensional 
interfaces. 

For analytical studies we employ simple power- 
counting arguments to establish the relevance or irrel- 
evance of the new terms at the equilibrium renormaliza- 
tion group fixed points. These will be combined with 
heuristic and scaling arguments, to relate the emerging 
critical behavior to equilibrium wetting and short-ranged 
non-equilibrium unbinding. 

For numerical studies, we consider one-dimensional 
discretizations of Eq. Q). As direct integrations of KPZ- 
like equations are known to be plagued with numerical 
instabilities resort to the exponential or Cole- 

Hopf transformation, n = exp(±/i), that leads to well- 
behaved, numerically tractable, Langevin equations with 
multiplicative noise j5i .IS] . In order to integrate these 
equations we eniploy a recently proposed efficient nu- 
merical scheme specifically designed to deal with 
stochastic equations with non-additive noise. More than 
just a useful technical trick, this transformation has an 
interesting physical motivation, as we discuss next. For 
negative values of a — Cc, the average interfacial height 
(h) (thickness of the liquid film) may be large but fi- 
nite, and the interface fluctuates around its average po- 
sition, occasionally touching the substrate. As the inter- 
face moves to infinity when a — > Oc, its average height 
grows (i.e. the liquid film completely wets the substrate) 
thereby suppressing contact (dry) sites. An appropriate 
order-parameter (OP) for the unbinding transition is the 
number of contact (dry) sites 0, , or equivalently the 
surface order parameter . This OP is finite and posi- 
tive when the interface is bound, and vanishes at the un- 
binding transition. The variable, (n) = (exp(— /i)), that 
vanishes exponentially far from the wall, is an adequate 
mathematical representation of such an OP (though not 
the only one). 

The main goal of our study is the description of the 
scaling behavior of the OP. (n) , is expected to obey sim- 
ple scaling near the critical point, for sufficiently large 
times, t, and large system-sizes L. Denoting 5a = |a— ad, 

{n{6a, t, L)) = L''^^^ {n{L^/'' 5a, L-H)), (5) 
while right at the transition {n{5a — 0,t)) ~ f-Pop/vz ^ 



l-6op g^j^(j therefore {n{5a, t = 00)) 5a^"'' , where 
the critical exponents were introduced following standard 
nomenclature. Analogously, for the interfacial height we 
can define {h) ~ 5a-'^>^ and {h{5a = 0,t)) ~ t'^'^/"'' - 
although in terms of /i a single universality class, with 
exponents related to the free KPZ is observed for 
both signs of A. Determining all of these critical expo- 
nents by the aforementioned techniques will allow us to 
assign the emerging critical behavior to specific univer- 
sality classes, providing a comprehensive classification of 
non-equilibrium unbinding transitions in the presence of 
long-ranged forces. 

Before proceeding to the presentation of our results, 
we notice that it is expected that the behavior for short- 
ranged interactions is recovered in the large-p limit of the 
long-ranged ones. Next, a brief review of the former is 
provided. 

III. BRIEF REVIEW OF NON-EQUILIBRIUM 
SHORT-RANGED UNBINDING 

The KPZ equation with exponential bounding poten- 
tials is 

dth = V^h + \{Vhf +a + b e'^^ + c e""'* + err], (6) 

with q > p > 0. The results for the four possible physical 
situations are: 

1. Repulsive wall and A < 

If A < (we set A = —1) the change of variables n = 
exp(— /i) transforms lO (with c = 0) into 

dtn = V^n — a n — b ti^^^ + narj. (7) 

This describes complete wetting transitions (along path 
1 in Fig. [2Ia)) characterized by (see 0) a dynamic expo- 
nent z = 3/2, identical to KPZ, = l/{2z - 2) = 1, and 
non-trivial exponents (3op and 6op that were determined 
by simulations. The exponents for h have been measured 
also and the transition was shown to be in the multi- 
plicative noise 1 (MNl) universality class: Pop ~ 1.78, 
Oop w 1.18, Oh w 0.33 and I3h « 1/2 in d = 1 ||. 

2. Repulsive wall and A > 

As for the A < case, it is more convenient |23 to use 
the transformation, n = exp(+/i) leading to 

dtn = V^Ti + art + bn^^^ + narj. (8) 

This equation describes the transition along the path 1 
in Fig. El^b). Numerical estimates for the associated 
universality class have been recently obtained from this 
(non-order-parameter) Langevin equation 20]. By mea- 
suring the order parameter m — (l/n) (that vanishes at 
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tive dynamics of the particle-like interface-surface con- 
tact points (i.e. points trapped in the potential well). 
Its leftmost border corresponds to the abrupt (discon- 
tinuous) binding of initially unbound interfaces. Again 
we refer the reader to 5] for a detailed discussion and 
to |23j for a review on generic phase-coexistence in non- 
equilibrium systems. 



Jf.. Attractive wall and A > 



FIG. 2: Phase diagrams for A < (a), A > (b). Paths la- 
beled 1 correspond to non-equilibrium complete wetting tran- 
sitions; 2, critical or first-order unbinding transitions (not 
studied); 3, first-order unbinding transitions; 4, unbinding 
transition in the directed percolation universality class. For 
A < and b < (attractive substrates) , two-phase coexis- 
tence is observed in the area delimited by the two lines. 



the transition), the following set of exponents was ob- 
tained: eop w 0.22, Pop « 0.32, different from MNl 
and 9h w 0.33, Ph ~ 0.5, z = 3/2 and = 1 in line with 
the corresponding exponents of the MNl class [M l20ll2ll |. 
This universality class is known as multiplicative noise 2 
(MN2). A detailed discussion of the differences between 
the MN2 and MNl universality classes, may be found in 
i 

Note that, apart from the signs, the difference between 
Eq. O and Eq. ||HJ) is in the leading power of n. It 
is possible, however, to summarize these two Langevin 
equations in 

dtu — V^n + aan + abn^+ narj, (9) 

with a = A/|A| and 7 = 1 — ap. Then a > 1 and a < 1 
correspond, respectively, to the MNl and MN2 universal- 
ity classes. In the first case the leading power for large 
values of n is the non-linear term while this role is taken 
by the linear term in the second case. The transition at 
the boundary 7 = 1 (p = 0) is obviously discontinuous, 
as both terms are linear and there is no saturating term. 

In MNl the order parameter is n, while in the MN2 
case, it is m = 1/n. In both cases a is the control pa- 
rameter. 



3. Attractive wall and A < 

For attractive walls, 6 < 0, a positive value of c is re- 
quired for stability, for any value of A. In systems with 
A < (see Fig. 1(a)), a new phenomenology including 
a broad coexistence region, and a dircctcd-percolation 
unbinding transition emerges 0,11 122, 23J. In the broad- 
coexistence region the stationary solution is either bound 
or unbound depending on the initial conditions 0, [3- 
Such a region is delimited on the right (where the bound 
phase loses stability) by a directed percolation transition, 
where the scaling properties are controlled by the effec- 



For A > (see Fig. 1(b)) a first-order transition sep- 
arates bound from unbound phases (akin to the equilib- 
rium discontinuous transition for attractive walls). No 
broad coexistence region, nor directed percolation tran- 
sition, exist in this case. 

IV. NON-EQUILIBRIUM LONG-RANGED 
UNBINDING: RESULTS 

We are now set to discuss the long-ranged non- 
equilibrium problem described by Eq. There is a 
singularity at /i = and thus only positive values of h 
are allowed (mimicking the impenetrability of the sub- 
strate). As before, if 6 > we take c = for simplicity. 
Proceeding as in the short-ranged non-equilibrium case, 
we perform the change of variables n = exp^ah), with 
a = A/|A|, in equation (0J, obtaining 

dfU = V^n -I- aan + ab- — - — — h ncrry, (10) 
|Q!log(n)|^+P 

where a term +acn/\a\og{n)\^'^'^ has to be added when 
6 < 0. As before, for positive A {a — 1), the order- 
parameter is m = 1/n, while for A < the order- 
parameter is n itself. Note also that as there is a sin- 
gularity at n = 1 (inherited from the singularity at h — 
in Eq. Q), for A > 0, where n diverges at the tran- 
sition, the initial condition is fixed at n{x) > 1 Va;, 
while for A < 0, where n vanishes at the transition, 
< n{x) < 1 Va;, is taken. The deterministic one-site 
terms of Eq. I|1U|) may be written as minus the derivative 
of a potential, U (n) , that is depicted in Fig. |21 

It is instructive to compare this model with the two 
universality classes reported for non-equilibrium short- 
ranged wetting, i.e. MNl and MN2. In fact, it is ex- 
pected that, in the limit of sufficiently large p, the power- 
law force yields the same dynamics as short-ranged (ex- 
ponential) forces. Thus, for A < and large p we an- 
ticipate MNl behavior while MN2 scaling should obtain 
when A > 0, in the same limit. 

A. Analytic results 

In an early work the KPZ non-linearity was argued 
to be irrelevant above the (mean-field) wetting tempera- 
ture 6 = &u) = 0, and an equilibrium (complete) wetting 
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FIG. 3: Effective potentials in terms of n obtained from tlie 
numerical integration of the interaction part of eq. Left 
and right panels correspond to, respectively, repulsive (6 — 1) 
and attractive (6 = —1) interactions. For both b we show 
results for negative and positive A, and different values of a 
corresponding to the bound and unbound phases, as well as 
at coexistence. Note that for A < 0(> 0) the unbound phase 
corresponds to a minimum of n at (oo). Transitions for at- 
tractive walls can be either first-order or continuous (directed- 
percolation) . 



transition was predicted to occur as a — > Oc, at constant 
b > bw, for any A ,24] (transitions along path 1 in Fig. 
121) . In the following we show that such a prediction is un- 
tenable and that the non-equilibrium term leads to new 
physics. 

Let us start by employing naive power counting argu- 
ments, based on equilibrium scaling, to decide whether A 
is a relevant or an irrelevant perturbation, at the mean- 
field fixed point and at the fluctuation one. In order to 
do that, we first fix A = in equation If 6 > 0, 
then the upper critical dimension depends only on the 
repulsive part of the potential and is dc{p) — 2p/(2 + p) 
[25. 26.] , as shown in the appendix. Now, from dimen- 
sional analysis [A] = Upon evaluating it at 
dc{p) one finds [A] = i^-2/(2+p) -vvhich, in terms of mo- 
menta, has a positive dimension for any value of p. There- 
fore, the KPZ non-linearity is relevant at the mean-field 
equilibrium wetting transition. 

Relevancy at the fluctuating regime fixed-point is 
proved using the known one-dimensional scaling dimen- 
sion of the field [h] ~ t^^^ ~ L^^^ at the fluctuation- 
dominated fixed point (see appendix). Then, it follows, 
[A] = implying that A is strongly relevant in 

any space-dimension. To be rigorous we would need to 
include perturbative corrections generated by the new 
non-linear term proportional to A, but even without com- 
puting these, one can say that it is very unlikely that such 
corrections reverse the strong lowest-order relevancy of A. 
The relevancy of A, is strongly supported by the results 
of numerical simulations of the corresponding Langevin 
equation as we will show next. 

As in one-dimensional equilibrium interfaces, where 
p — 2 separates the mean-field and the fluctuation- 
dominated regimes, it is easy to argue that in non- 
equilibrium the two regimes are separated by p = 1. 



From Eq. Q in the absence of noise, the mean-field ve- 
locity exponent at the critical point, given by A((Vft.)^) -I- 
flc = 0, is obtained by integrating dth ^ h^P^^, and 
found to he 9h = l/{p + 2). On the other hand, when 
noise (fluctuations) is included, the (one-dimensional) 
free KPZ equation has a roughening exponent of 1/3 
and, therefore, a velocity proportional to t^^^ "oj. Which 
of these contributions dominates? Clearly, if p < 1 the 
wall-induced velocity is larger and fluctuations give only 
a higher order correction (i.e. they are irrelevant). By 
contrast, if p > 1 the effective repulsion generated by 
the wall (through suppression of the intrinsic interfacial 
roughness) controls the scaling. Thus, in non-equilibrium 
long-ranged wetting, p = 1 separates the mean-field from 
the fluctuation-dominated regimes. 

Transient effects, that are significant before the non- 
equilibrium interface develops its full (asymptotic) time- 
dependent roughness, may prevent the KPZ exponent 
= 1/3 from being observed, leading to an effective ex- 
ponent, 9eff < 1/3. Furthermore, at short times, the 
interface is expected to grow with an Edwards- Wilkinson 
exponent, 9 = 1/4, and therefore 9eff, increases progres- 
sively from 1/4 to its asymptotic KPZ value, 1/3, in the 
long time regime. Comparing these values with the wall 
induced velocity exponent l/{p + 2), we anticipate that 
for potentials with 1 < p < 2 severe transient effects will 
occur before the fluctuation-dominated scaling sets in. 
By contrast, for p > 2 fluctuations dominate from the 
early stages of interfacial growth. 



B. Numerical Results 

In order to avoid numerical instabilities, typical of KPZ 
direct numerical integration schemes [15i |. we chose to 
study the associated multiplicative noise Eq. H1U|I ob- 
tained after performing a Cole-Hopf transformation. To 
solve Eq. (|10|) efficiently we have used a recently pro- 
posed split-step scheme for the integration of Langevin 
equations with non-additive noise UT'I. In this scheme, 
the equation under consideration is discretized in space 
and time and separated into two contributions: (i) the 
first includes deterministic terms only and is integrated 
at each time-step using a standard integration scheme: 
Eulcr, Rungc Kutta, etc |23| (here we have chosen a sim- 
ple Euler algorithm) (ii) the output of the first step is 
used as input to integrate (along the same discrete time- 
step) the second part which includes the noise and, op- 
tionally, linear deterministic terms. This is done by sam- 
pling the probability-distribution, i.e. the solution of the 
Fokker-Planck equation associated with this part of the 
dynamics. In the case under study (noise proportional 
to the field), the second step can be carried out exactly. 
At each site, one has to sample a log-normal distribution, 
i.e., the solution of the Fokker-Planck equation associated 
with dfU = aan + anr] (for more details see ||2^ and jl7i|). 
The two-step algorithm for Eq. (|10|) is then implemented 



6 



as follows. At each site n = n{x,t), we compute 

ab n „ 



ni{x,t) ~ n + dt 



l+p ' ^ discr 



n{x, t) 



11) 



(Q!log(n)) 

where the discretized Laplacian is defined by 

2 f n{x + Ax,t) + n{x — Ax,t) — 2n{x,t) 



"^MscMx, t) 



with Ax the space-mesh, and 



(12) 



n{x^t + At) — ni[x,t) exp ( aaAi + (T 77 \/At J (13) 



where 77 is a random variable extracted from a Normal 
distribution with zero-mean and unit variance. Note that 
the linear deterministic term can be included in either the 
first or the second step, or partially incorporated in both 
of them. For systems with 6 < 0, the stabilizing term, 
proportional to c has be to included. 

We set (7 = 1, Aa; = l/VO.l, and the time- mesh 
At = 0.1 (note that in this scheme At can be taken 
larger than in the usual integration algorithms In 
some simulations we used different values of 6, which by 
default was set to 6 = ±1. We take as initial condition 
n{x,t = 0) = 3 if A > (recall that n e]l,cx3[) and 
n = 0.5 if A < (n e [0, 1[). Then, the dynamics is iter- 
ated by employing the two-step integration algorithm at 
each site and using parallel updating. 

The numerical procedure is as follows. In order 
to determine the critical point for any set of parame- 
ters we take the system-size as large as possible and 
look for the separatrix between upward-bending and 
downward-bending curves in the order-parameter (either 
noT m = 1/n depending on the case) versus t in a double- 
logarithmic plot. The asymptotic value of this slope gives 
an estimation of 9op- Also, for the same parameters, 
(h) grows as a power-law with an exponent 0h (bending 
downward and upward in the bound and the unbound 
phases, respectively). Generally the order parameter is 
more sensitive to control-parameter variations, providing 
the most reliable way of determining the critical point. 
For completeness, and in order to check the validity of 
analytical approximations, we measure the global inter- 
face width, W, at the transition, which is expected to 
grow with the KPZ exponent (3w = 1/3, in the regime 
where it is asymptotically free. 

Once the critical point is determined accurately we 
compute Po p and j3h by measuring the stationary values 
of the order parameter and of {h) at different distances 
from it. A complementary approach is based on finite- 
size scaling analysis: the values of the order parameter 
and of {h), at saturation, are measured for a fixed value 
of a as a function of system-size. At the critical point 
these values scale with exponents Pop/v and Phi v., re- 
spectively. In addition, the scaling of the saturation times 
for different system-sizes allows to determine the dynam- 
ical exponent z. This standard finite-size scaling analysis 



is not always possible (see below), and in such cases z is 
measured through spreading-like experiments. Finally, 
an alternative to spreading consists in measuring the dis- 
tribution of gaps between contact points at a given time. 
For small gaps this function decays with an exponent 
zOq p giving yet another estimate of z 29] . 

The correlation length critical exponent v is obtained 
by measuring the location of the effective critical point, 
i.e. the value of a for which the order-parameter falls 
below a fixed threshold, say 10~^, as a function of system- 
size; ac^eff{L) ^ L^^/'^ . This exponent may also be 
determined indirectly by employing scaling relations and 
using the value of Pop Iv from finite-size scaling analysis 
and the value of Pop obtained from direct measurements. 
The results of these measurements, in conjunction with 
the scaling laws, provide an over-complete estimation of 
the set of critical exponents, that was also used to verify 
scaling relations. 

Before discussing the differences between the various 
universality classes and regimes (i.e. different values of 
A, p and b) we first give an overview of the common 
features of all simulations. 

1. Once the KPZ equation parameters {D,X,a) are 
fixed, the location of the critical point is univer- 
sal, meaning that it does not depend on the details 
of the substrate, i.e. on the values of 6, c, and p. 
The critical point is determined by the value of a 
where the free-KPZ interface changes the sign of 
its velocity, from positive, i.e. diverging to an un- 
bound state, to negative, becoming bound at the 
wall: Oc -I- A((V/i)^) = 0. As we consider two differ- 
ent values of A, +1 and —1, there are two critical 
points: a^X = ±1) ~ tO. 143668(3). 

2. At the critical point, the asymptotically unbound 
interface is a free KPZ one, and thus z = 3/2 and 
Pw = 1/3. These values were consistently checked 
in all simulations (see Fig. EIa),(b)). 

3. A simple argument, originally given in 'l6l, predicts 
V = 1 for all bounded KPZ interfaces. This predic- 
tion was confirmed in all of our cases (see inset of 
Fig. Ha)). 



1. Repulsive walls and A > 0. 

We have to distinguish two regimes, depending on the 
range of the attractive substrate, i.e. the value of p. 

Mean-field regime. The theoretical discussion indicates 
that for p < 1, and any sign of A, a mean-field regime 
controlled by the exponents Oh = l/{p + 2) and Ph — 
3/(2p + 4) is obtained. By changing variables in a 
naive way a stretched exponential behavior for the order- 
parameter is predicted. Figures [S] and ^| illustrate the 
confirmation of these predictions (both for positive and 
negative A). 
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FIG. 6: Log-log plot of the time evolution at ac of (h) (dashed, 
black curve), and — log{nop) (solid, red curve), in the fluc- 
tuation regime (p = 2) and for A > 0. From the slopes of 
the straight-line fits one finds 9h = 0.32(2) (upper-curve) and 
9op = 0.228(6) (lower-curve). Upper-inset: finite-size scaling 
of (h) yielding jSh/v = 0.46(2). From the lower-inset one ob- 
tains Pop/v = 0.34(2). These exponents agree with those of 
the MN2 universality class. 



FIG. 4: Features common to all simulations, (a) Roughness 
vs t gives (iw ~ 0.33(1); (inset) —ac{L) vs lO'^L"^ falls on a 
straight-line that yields ~ 1 (data for A = — 1 and p = 2). 
(b) Saturation time vs system-size leads to z = 1.48(4) (data 
for A = 1 and p = 2) . 
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FIG. 5: Log-log plot of the time evolution at of (h) (upper, 
red curves), — ln{nop) (middle, green curves), and the width 
w (lower, black curves), in the mean-field regime p = 0.5, 
for A < (main) and A > (inset). Irrespective of the sign 
of A, (h) and the roughness may be fitted to a power-law 
with the predicted exponents 9h ~ l/(p -I- 2) and f3w = 1/3, 
respectively. — In(nop) falls on a straight-line in a double log- 
arithmic plot, confirming the stretched-exponential behavior 
of the order parameter. 



Fluctuation regime: Multiplicative Noise 2. A strong- 
fluctuation regime is predicted for systems with p > 1 
but, as argued above, severe transient effects are expected 
for 2 > p > 1. We start with the analysis of the, a priori, 
simpler p > 2 sub-regime and offer simulation results for 
p — 2,2.5,3,4,7. In all cases the order-parameter was 
found to decay at criticality with an exponent 9op ~ 
0.229 while the average height diverges with 6h « 1/3 
(see Fig. |S1 data shown for p = 2). A standard flnite-size 
scaling analysis can be performed (see Fig. EJ, yielding 
Pop/i, = 0.34(2) and Ph/i^ = 0.46(2). 

These results, together with the previously reported 
general ones, unambiguously place the fluctuation regime 
for repulsive walls with positive A into the MN2 univer- 
sality class. 

For systems with 1 < p < 2, where strong transients 
are expected, after fixing b — 1 and running simulations 
up to t = 10^, continuously varying power-law exponents 
are found (see Fig. [T)). We note, however, that these fits 
give effective rather than asymptotic exponents. In fact, 
the change in the effective exponents from mean-field 
(wall-controlled) to the fluctuation (intrinsic-interface) 
regime is expected to occur at shorter times when the 
effect of the substrate is less pronounced, implying that 
reducing the value of b decreases the crossover time. This 
was confirmed by simulating systems with b = 0.1 and 
b — 0.05 and observing a monotonic decrease of the effec- 
tive exponents that converge to the expected asymptotic 
value 6'op ~ 0.228, 6*/, « 1/3 (see inset (a) of Fig. d) 
in line with the hypothesis that the transition belongs to 
the MN2 universality class. 

In order to check that p = 2 is the boundary between 
the strong and weak transient sub-regimes, we have plot- 
ted in Fig. \7\ inset (b), the average order-parameter for 
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FIG. 7: Main: Order-parameter multiplied by the expected 
power- law f"'^^*. For systems with 6=1, long transients that 
depend on p are observed for systems with 1 < p < 2. Inset 
(a): The crossover times are reduced as h decreases; compare 
the upper, red curve for h = 0.05 with the lower, black one 
for h — 1 {p = 2). Inset (b): Order-parameter at t = 10^ vs 
p {b = 1). p — 2 marks the boundary of the strong transient 
region as illustrated by the different behaviors observed above 
and below p = 2. 



systems with the same initial condition, at time t = 10^, 
and different values of p. This is a non-stationary value 
of the OP that is strongly affected by transients. It is 
clear from the figure that the behavior of the order pa- 
rameter changes qualitatively at p = 2 corroborating the 
result that this value of p marks the boundary between 
the sub-regimes with and without severe transients. 

2. Repulsive walls and A < 0. 

Mean-field regime. In parallel with the positive A case, 
the results of Figs. Isl and [TUl show that the theoretically 
predicted mean-field regime is clearly observed for sys- 
tems with p < 1. 

Fluctuation regime: Multiplicative Noise 1. Again we 
have to distinguish two sub-regimes, with and without 
severe transients, depending on whether p is larger than 
or smaller than 2. Simulations in the weak transient 
regime were performed for p = 2 and 3. In both of these 
systems the order-parameter decays at criticality with 
an exponent Oqp — 1.19(1) while the average height di- 
verges with 9h — 0.33(1) (see Fig. |lta),(b) data shown 
for p — 2). As was first pointed out in '2§\ , finite-size scal- 
ing measurements are non trivial in this case due to the 
presence of two different characteristic times. Namely, 
the correlation length reaches the size of the system at 
times ~ L^, whilst the interface typically detaches from 
the wall at times ~ L^/^". As the latter grows with a 
larger exponent for MNl, the interface detaches from the 
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FIG. 8: Results for A < and p = 2. (a) Log-log plot of 
the time decay of the order-parameter. The straight-line is a 
guide to the eye, and has a slope 9op ~ 1.19(1). (b) Log-log 
plot of the average distance from the wall vs time, leading to 
9h ~ 0.33(1). (c) The scaling of the saturation value of the 
order-parameter yields /3op ~ 1.76(3). (d) Log-log plot of 
the gap (between contact points) distribution function. The 
initial slope gives z6op ~ 1.75(10). 



wall before it reaches the saturation regime for finite sam- 
ples, rendering the evaluation of Pop/^ and z through 
standard finite-size scaling methods problematic. An es- 
timation of (3op is possible by taking a large system-size, 
L = 2^^, and measuring the order-parameter stationary- 
state value upon approaching Oc. We find Pop = 1.76(3) 
and (3h — 0.51(3)(see Fig. EIc)). z is accessible through 
spreading experiments from an initial condition with only 
one active (pinned) site. The measurement of the mean- 
square deviation from the origin R^{t) « i^^^, gives 
z = 1.52(5) (not shown). Alternatively, one can in- 
vestigate the gap distribution function of the distances 
between neighboring contact points at a given time or, 
equivalently, the average size of inactive islands in the n- 
language [23 ■ For small gaps this function decays with 
an exponent zOqp, and we find zOqp — 1.75(10) which 
leads to a value of z compatible with 3/2 (see Fig. |HJd)). 

These results, together with the general ones, place 
unambiguously this fluctuation regime for repulsive walls 
with negative A into the MNl universality class. 

Again, for systems with 1 < p < 2, different effec- 
tive exponents are obtained at a fixed maximum time for 
different values of b (unity and smaller), confirming the 
existence of strong transients. Upon decreasing 6, the in- 
fluence of the wall is reduced and a behavior compatible 
with the MNl class is observed: Oqp ~ 1.19, 9h ~ 1/3, 
Pop ~ 1.76 and Ph ~ 0.5 (figure not shown). 



3. Attractive wall and A < 0. 

The phase-diagram, depicted on the left panel of Fig. 
is similar to that found for short-ranged interactions 
. For a fixed b, by varying a one of two transitions 
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may occur depending upon the initial interfacial state. 
Initially unbound interfaces experience an unbinding- 
binding transition at Oc where the free-interface velocity 
inverts its sign (in full analogy with the previous case; see 
path 3 in Fig. On the other hand, initially bound in- 
terfaces unbind at a different non-trivial value of a, noted 
a* > Oc, inside the free- interface unbound phase (path 4 
in Fig. |2Ia)). This transition is analogous to the one ob- 
served for short-ranged forces, and is expected to be con- 
trolled by the unbinding of interface-sites trapped in the 
potential minimum. Bound sites (located around the po- 
tential minimum) are identified with particles; unbound 
sites are described by holes. The effective particle 
namics is very similar to that of the contact-process |23 
(a well-studied model known to be in the directed perco- 
lation class): an occupied site can become empty when 
a point is detached, and can induce also the binding of 
a neighboring site. Furthermore, empty sites cannot be- 
come spontaneously occupied in the absence of occupied 
(bound) neighboring sites. Indeed, as soon as the inter- 
face is locally out of the potential well, it is pulled away 
from it. This corresponds to the absorbing state charac- 
teristic of the directed percolation class. Note that the 
statistics of the average number of such pseudo-particles 
is completely analogous to that of (exp(— /i)). 

Before the depinning transition, typical triangular 
structures are observed, consisting of pinned sites (lying 
in the potential well), and depinned sites being pulled 
from the substrate. This triangular shapes (pyramidal 
in two-dimensions) are similar to those in the analo- 
gous short-range case, and are reminiscent of pyrami- 
dal mounds obtained in the non-equilibrium growing of 
some interfaces, as for instance, in the so called Stranski- 
Krastanov effect . 

Our numerical results show that this transition is con- 
trolled, as in the short-ranged case, by directed perco- 
lation critical exponents (see Fig. inj. In particular, we 
have determined Pop/i" = 0.26(2) and Oop = 0.161(2), 
in excellent agreement with the one-dimensional directed 
percolation values. 

Let us remark that, as the bound sites remain inside 
the potential well, and the dynamics controlling their fi- 
nal "escape" is likely to be insensitive to the exponen- 
tial or power-law tails of the potential at large values of 
h, the parallel between this behavior and the directed- 
percolation transition for short-ranged forces is to be ex- 
pected. Interestingly, Ginelli et al. investigated a lattice 
model of a generalized contact process with long-ranged 
interactions between the edges of low-density segments 
and found a transition in the directed percolation univer- 
sality class for forces that decay sufficiently slowly, and a 
first-order transition otherwise |3 1| . Clearly, in terms of 
h this translates into a long-ranged interaction between 
the vertices forming the triangle bases, and it is reason- 
able to assume that, in turn, an effective long-ranged 
attraction between the substrate and the interface must 
be obtained. In the light of these results it is reasonable 
to assume, that both short- and long-ranged interactions 
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FIG. 9: Time decay of the order-parameter at the critical 
point a = 0.38448. From the slope in the log-log plot we find 
9 — 0.161(2). Inset: finite size scaling of the order parameter. 
From the slope in the log-log plot we estimate P/v = 0.26(2), 
in agreement with directed percolation values. 



in similar models will be characterized by the same be- 
havior below bw 

In the region between Cc and a* one observes generic 
(broad) phase coexistence: the stationary solution is ei- 
ther bound or unbound depending on the initial condi- 
tion. Within this region, the bound phase is character- 
ized by some bound sites trapped in the potential mini- 
mum, and pseudo- unbound regions separating them J]. 
In full analogy with short-ranged forces, close to the un- 
binding transition a < a* initially bound interfaces are 
stable owing to a mechanism that eliminates local fluc- 
tuations into the unbound phase: once formed, islands 
of the unbound phase rapidly transform into triangular 
mounds of fixed slope, which subsequently shrink from 
the edges. 



4-. Attractive wall and A > 0. 

When A > 0, the situation is rather similar to the one 
for equilibrium and for non-equilibrium (A > 0) short- 
ranged systems. At the critical value ac where the free 
interface inverts its velocity sign, there is a discontinuous 
unbinding-binding transition (path 3 in Fig. [^J. This 
value does not depend on the value of p nor on 6 or c 

M 

The multicritical point. Finally, for either sign of A, path 
2 in Fig. [21 corresponds to a multi-critical point analog to 
an equilibrium critical wetting transition when the crit- 
ical point is approached at coexistence. Most likely, its 
location will not coincide with its mean-field value 6 = 0, 
but exhibits some renormalization shift. The analysis of 
this multi-critical point will be considered elsewhere. 
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FIG. 10: Time growth exponents Oh in the mean-field like 
regime for (/i) at the critical point, as results from the discrete 
interfacial model. Blue circles(red squares) stand for A > 
0(< 0) data points, and the solid line is the predicted curve 
l/(p + 2). 



C. Discrete Model 

As a final check of universality issues, we simulated a 
discrete interfacial model, known to belong to the KPZ 
class, in the presence of a long-ranged substrate. The 
model is the same as that studied in the context of short- 
ranged wetting in |2]| . Even if plagued with long tran- 
sient effects (much larger than in the short-ranged case) 
all of the previously reported phase diagrams and univer- 
sality classes seem to be confirmed for the different types 
of walls (i.e. values of h and p) and signs of the non- 
linearity. Generally, the discrete model provides slightly 
better results for the height variable as compared with 
the continuum model, and worse for the order-parameter. 

Figure ^| displays the time growth of the mean sep- 
aration {h) in the mean-field like regime (p < 1), for 
both positive and negative A. Additionally, the ratio 
(3op/v = 0.251(2) and 9op = 0.156(2) were obtained 
for the directed percolation transition, which compares 
favorably with the accepted estimates 0.25208(5) and 
0.1595 [33. 

V. DISCUSSION AND CONCLUSIONS 

We have studied the unbinding of KPZ interfaces in 
the presence of limiting substrates, interacting via long- 
ranged potentials. This is the simplest model for interfa- 
cial effective descriptions of wetting and in general, un- 
binding transitions, of systems interacting through van 
der Waals forces under non-equilibrium conditions. 

We have presented the results of systematic analyti- 
cal and numerical studies of one-dimensional KPZ-like 
interfaces in the presence of long-ranged forces Eq. Q , 
supporting the following conclusions: 

i) Repulsive interactions drive a non-equilibrium com- 



plete wetting transition for either sign of A. This tran- 
sition belongs to different universality classes depending 
on the strength of the repulsion, i.e. on the value of p 
in Eq. Q and on the sign of A. For p < 1 a mean-field 
like regime is observed in both cases, while for p > \ 
the fluctuation regime obtains and the transition is in 
the multiplicative noise 1 (MNl) class for A < and in 
the multiplicative noise 2 (MN2) for A > 0. Systems 
in the fluctuation regime, exhibit severe crossover effects 
for bounding potentials with \ < p < 2. This should 
be contrasted with the behavior of equilibrium systems 
where the value of p that separates the mean field from 
the fluctuation regimes was found to be p = 2. More im- 
portantly, in non-equilibrium systems the symmetry of 
the equilibrium wetting and drying transitions is broken 
and the fluctuation regime of the corresponding equilib- 
rium wetting transitions is split into two different non- 
equilibrium universality classes, MNl and MN2 respec- 
tively. Our results are collected in table II II 

ii) For attractive walls, i.e. below the critical wetting 
temperature, phase-diagrams analogous to those of sys- 
tems with short-ranged forces have been found: generic 
phase-coexistence over a finite area limited by directed 
percolation and first-order boundary lines for A < 0, 
and a first-order phase transition from an unbound to 
a bound interface for A > 0. This transition should not 
be called "wetting" as the interface detaches below the 
wetting transition temperature. 

The unbinding transition at the critical wetting point 
(which in the language of this paper corresponds to a 
multicritical point) requires a higher degree of fine-tuning 
and is therefore expected to be more difficult to observe in 
experimental situations. Its study is also more laborious 
and is deferred to future work. 

For more realistic two-dimensional interfaces, corre- 
sponding to three-dimensional bulk systems, the situa- 
tion is expected to be very similar: all universality classes 
(mean- field, multiplicative noise 1, multiplicative noise 2, 
and directed percolation) are expected to be substituted 
by their two-dimensional counterparts, with analogous 
phase diagrams and overall phenomenology. 

We hope that the results described in this paper will 
help to motivate an experimental study of wetting and 
unbinding transitions under non-equilibrium situations. 
In these systems one expects to find the rich phenomenol- 
ogy described here, and they can be used to test some 
of our quantitative predictions, concerning the values of 
the exponents and the existence of various universality 
classes. Liquid-crystals "g], molecular-beam epitaxial sys- 
tems, as GaAS 7] claimed to grow following KPZ scaling, 
or materials exhibiting Stranski-Krastanov instabilities 
8] , appear to be good candidates that are at least worth 
investigating in this context. Indeed, it is rather excit- 
ing to think that non-equilibrium complete wetting ex- 
ponents are measurable. This would be a way of measur- 
ing the multiplicative noise critical exponents, and brings 
new hope of measuring directed percolation exponents in 
real systems [s^- 
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FIG. 11; Summary of the critical exponents in the mean-field (p < 1) and the fluctuation (p > 1) regimes for non-equilibrium, 
complete wetting transitions with long-ranged forces. To facilitate the comparison, the exponents for the MNl and MN2 
universality classes are also included (p = oo). *, exponent from finite-size analysis or scaling relations; ^, estimated value from 
short simulations; n.a., not applicable. 
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APPENDIX: BRIEF REVIEW OF EQUILIBRIUM 
WETTING 

The action associated with Eq. |^ j2^ [s^ (setting 
c = 0) is 

S{h, h)^ J d'^xdt _ h ^Q^h - V^/i -a- bh-P-^] I , 

(A.l) 

where h denotes, as usual, the response field [26L l35l |. If 
one assumes first that the interaction term is the dom- 
inant one, from naive dimensional analysis, imposing b 
to be dimensionless at the upper critical dimension, and 
equating the dimensions of the time-derivative and the 
potential terms, one obtains [/i]mf = L'^/(-P+'^) and con- 
sequently, within mean- field, 6^ = l/{p + 2) since time 
scales naively as L^. The exponent values Ph = l/(p+ 1) 
and p = {p + 2)/{2p + 2) are then obtained by matching 
[a] = [h]Jfp and by identifying L with the characteristic 
correlation length, respectively. 

On the other hand, when fluctuations (i.e. the noise 
term) dominate, we require the noise amplitude to be 
dimensionless at the upper critical dimension, which 
leads to [h]FL = L^^+'^^/^, and therefore [hh] = L"'', 
[hjpL — L^^^'^'f^^. From this, proceeding as before, 
Oh = {2-d)/A,v = 2/{d + 2), and (3h = {2-d)/{d + 2). 
These results (which may be obtained u sing a number 
of different procedures [E S IH IM El Hj) are exact 
as long as h and h do not have anomalous dimensions, 
which has been shown to be the case js^ . 

The upper critical dimension is defined by [h]MF = 
[h\FL, which yields dc{p) — 2p/{2+p). Note that for d > 



dc{p) the critical exponents depend on the details of the 
interaction (i.e. on p) whilst for d < dc{p), they depend 
only on d. In particular, in one dimensional systems, 
p ^ 2 marks the transition between a mean-field regime 
and a fluctuation regime: 

i) If p < 2 mean field theory is valid, and consequently 
9h = 1/ip + 2), I3h = 1/ip +l),z = 2, and p = 1/2. 

ii) For p > 2, the substrate interaction decays fast 
enough for the fluctuations to take over and the expo- 
nents become p-independent: 9h = 1/4, /3;i = 1/3, z = 2, 
and p = 2/3. 

Note that at the limiting value p ~ 2 the exponents 
change continuously from the mean field to the fluc- 
tuation regime. It is also remarkable that the fluctu- 
ation regime exponents coincide with those of short- 
ranged equilibrium wetting (characterized by exponential 
bounding potentials Q). 

Until now we have considered the scaling properties of 
(/i), but as was mentioned earlier the number of dry sites 
or contact points between the interface and the substrate, 
measured by (exp(— /i)), is known to exhibit interesting 
scaling behavior in wetting problems |l9j| . 

i) For p < 2 simple mean field scaling holds, and the 
/i-distribution is a Gaussian detaching from the wall at 
a speed controlled by its mean value. As the interface is 
well described by its average position, it is expected that 

(e-'^) ^ e-W ^ 6-4*^''"^^', (A.2) 

yielding a stretched exponential decay. 

ii) For p > 2, {a + bexp{—h)) = holds in the station- 
ary state, and therefore (n) oc a: using simple scaling, 
[exp(-/i)] = [a] = [dth] ^ t^/^-\ giving (exp(-;i)) ~ 

xhis result can be derived in a number of ways, 
inclu ding explicit calculations for discrete models in this 
class p39| . and remains valid for long-ranged potentials 
in the fluctuation regime. Note the difference between 
this fluctuation-induced power-law behavior and the pre- 
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viously reported stretched exponential behavior in the 
mean-field regime. 

For attractive walls, b < 0, a positive value of c is 
required to ensure the impenetrability of the substrate 
(see Fig. In this case it is easy to argue that the 
interface jumps discontinuously from a bound state (for 



a < 0), localized at the minimum of V{h), to an unbound 
state (for a > 0) through a first-order phase transition. 
Clearly, in terms of the contact points, (exp(— /i)), the 
transition is also discontinuous. 
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